Cascade approach to current fluctuations in a chaotic cavity 
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We propose a simple semiclassical method for calculating higher-order cumulants of current 
in multichannel mesoscopic conductors. To demonstrate its efficiency, we calculate the third and 
fourth cumulants of current for a chaotic cavity with multichannel leads of arbitrary transparency 
and compare the results with ensemble-averaged quantum-mechanical quantities. We also explain 
the discrepancy between the quantum-mechanical results and previous semiclassical calculations. 



I. INTRODUCTION 

In the last decade, there has been a large .interest in 
current correlations in mesoscopic conductors E Recently, 
also higher cumulants of current have received a aigjaif- 
icant attention of theorists. In a series of papersoo a 
scattering approach to the distribution of charge trans- 
mitted through an arbitrary multi-terminal, multi-mode 
mesoscopic conductor, i.e. the so-called full counting 
statistics has been developed. As was shown in Ref.B, 
the ensemble-averaged cumulants of arbitrary order can 
be calculated for any two-terminal conductor where the 
distribution of the transmission eigenvalues as known, 
e.g. for diffusive wiresR. chaotic cavities,Bu, double- 
barrier tunnel junctionsEHI, or combinations of different 
conductors . 

Recently, Nazarovtil presented a method for calculat- 
ing the full counting statistics of charge transfer in con- 
ductors with a large number of quantum channels based 
on equations for the semiclassical Keldysh Green's func- 
tions. Subsequently, this method was extendedtj to mul- 
titerminal systems. Also, other approaches_to higher cu- 
mulants, such as the nonlinear sigma modelij for diffusive 
wires, have been proposed. _ __ _ 

Common to all approacheso □'EtE^I is that they are 
based on a quantum mechanical formulation. To ob- 
tain the cumulants for semiclassical systems, i.e. systems 
much larger than the Fermi wavelength, an ensemble av- 
erage is performed and the number of transport modes is 
set to infinity, i.e. single-mode weak-localization-like cor- 
rections are neglected. Therefore it is of interest to have 
a completely semiclassical theory for the higher cumu- 
lants, which does not involve any quantum-mechanical 
quantities. _ 

A step in this direction was made by de Jongcj, who 
calculated the distribution of charge transmitted through 
a double-barrier tunnel junction by applying a master 
equation to the transport in each completely indepen- 
dent transverse quantum channel. The results were in 
agreement with the quantum-mechanical theory in the 
limit of large channel number. 

An attempt to construct a fully semiclassical theory 
of higher cumulants of current in a chaotic caidty was 
made by Blanter, Schomerus, and BeqeakkerliJ based 
on the principle of minimal correlations .cJ According to 



this principle, the fluctuations of the semiclassical distri- 
bution function of electrons in the cavity and the fluc- 
tuations of outgoing currents are related only through 
the condition of electron-number conservation, which is 
equivalent to the dephasing-voltage-probe approaclJ13 in 
quantum mechanics. However an attempt to extend the 
minimal-correlation approach to the fourth cumulant has 
led to a discrepancy with quantum-iwechanical results,c£l 
which highly surprised the authors JIj 

Meanwhile the correlations imposed by the particle- 
number conservation are not the only possible ones 
in semiclassics. Quite recently the semiclassical 
Boltzmann-Langesan approaclJIa has been extended to 
higher cumulants. t£l This extension takes into account the 
effect of lower cumulants on higher cumulants through 
the fluctuations of the distribution function and there- 
fore it was termed cascade approach. Its equivalence 
with quantum-mechanical resultsu has been proven for 
diffusive metallic conductors. In this paper, we show 
that the cascade approach is not restricted to diffu- 
sive metals or to conductors where the scattering is de- 
scribed by a collision integral, but it may be also ap- 
plied to other systems that allow a semiclassical descrip- 
tion, e.g. to chaotic cavities. To this end, we semiclassi- 
cally calculate the third and fourth cumulants of current 
in a chaotic cavity taking into account cascade correla- 
tions and show that these values coincide with ensemble- 
averaged quantum- mechanical results. 

The paper is organized as follows. In Section |0| we 
describe the model of chaotic cavity to be considered. 
The minimal-correlation results for the second cumulant 



of current are presented in Section [II. In Section IV we 
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FIG. 1. A chaotic cavity with imperfect leads 
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calculate the third and fourth cumulants of current by 
means of the cascade approach. In Section [v], the 
third cumulant of current is calculated by means of the 
quantum-mechanical circuit theory and its equivalence 
with the cascade results is shown. Section [Vl] presents a 
conclusion where the results are summarized. 



II. THE MODEL 

Consider a chaotic cavity with two contacts of arbi- 
trary transparency. The left contact has Nl 3> 1 chan- 
nels and transparency Tl, and the right contact has N R 
channels and a transparency T R . The conductances of 
the leads G L = (e 2 /nh)N L T L and G R = {e 2 /irK)N R Y R 
are also assumed to be much larger than e 2 /h, and the 
total conductance pf, the system is that of two resistors 
connected in seriesEHI 



G = 



GlGr 
Gl + G} 



Because of strong chaotic scattering in the cavity, the 
electrons entering the cavity lose the memory of their 
phase on the time scale of the order of the time of flight 
through the cavity yet retain their energy. Therefore de- 
spite the quantum nature of the contacts, the cavity is 
a semiclassical object in the sense that the electrons in- 
side it may be described by a semiclassical distribution 
function that depends only on the electron energy. Its 
average value is given by an expressions 



m = 



G L f L (e) + G R f(e) 
Gl + Gr 



(1) 



where and f R (e) are the distribution functions in 

the left and right electrodes. 



III. THE PRINCIPLE OF MINIMAL 
CORRELATIONS 

If the distribution function in the cavity were not al- 
lowed to fluctuate, the cavity could be considered just 
as a reservoir with non-Fermian distribution of electrons. 
The contacts would be independent generators of cur- 
rent noise and the cumulants of the corresponding extra- 
neous noise currents could be obtained by differentiat- 
ing the cepjesponding quantum-mechanical characteristic 
functionsc 2 ] for the charge transmitted in time t 



XL,n{\t) = cxp 



fN 



L,R 



2nh 



cfeln{l 



+iW( £ )[i- f L A?)W x -i; 

+ T L ,Rf L Ms)[l~f(e)}(e~ lX -l)} 



(2) 



with respect to the parameter A the corresponding num- 
ber of times. 

In what follows, we will be interested only in the 
Fourier transforms of the current cumulants in the low- 
frequency limit and it will be implied that all the sub- 
sequent equations contain only low-frequency Fourier 
transforms of the corresponding quantities. Equation (j^) 
leads to the following expressions for the cumulants of 
the noise current generated by the contacts: 



M,r))= / de((I2 iR )) e , 



(3) 



where 



«^U»e = G L . R [!l.r(1 - f) + /(l - fL.R) 

- ^l,r(/l,r - f) ], 



(4) 



= eG L<R (f - f L<R ){l - 3T LtR [f(l - f LtR ) 



- /)] + 2r£ Jf - f L ) 2 }, 



(5) 



and 



((li R ))e = e 2 G L , R {f L , R {l - f) + /(l - f L . R ) 



+T L ,R(12f 2 f + 12f 2 f L , R - 12/1 / 2 



-VIr 



7f 2 + 2ff LtR ) 



+ur 2 L jf L> R - f) 2 [f L ,R.a - /) + /(i - fL, R )] 



6T 



l,r(/l,r - /) 4 }- 



(6) 



Since we are interested here only in low-frequency fluc- 
tuations, the pile-up of electrons in the cavity is forbid- 
den. On the other hand, the noise currents II and Ir are 
absolutely independent, which would apparently result in 
a violation of the current-conservation law if these were 
the only contributions to the current noise. To ensure 
the current conservation at low frequencies, one has to 
take into account fluctuations of the distribution func- 
tion Sf(e) in the cavityH Now the fluctuations of the 
current outgoing from the cavity to the left and right 
electrodes assume a form of Langevin equations, where 
II and Ir play the role of extraneous sources. 



6Il,r = Il,r + ~G l , r I ,l:dfc) 



(7) 



Extracting Sf from the condition of current conservation 
SIl + SIr = 0, 
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one obtains 



Sh 



GrIl — GlIr 
Gl + Gr 



(8) 



By squaring this equation and using the independence of 
II and Ir, one easily obtains that the second cumulant 
of the measurable current is 



_ <%<<%))+<%((%)) 

[[ L>> (g l + c R r 



(9) 



In the zero-temperature limit it gives 

((/£)) = e/ [GlGr(Gl + G L ) + G|(l - r fl ) 

+ G*(l - r L )] /(G L + G*) 3 , (10) 

where I is the average current flowing through the cav- 
ity. In the high-transparency limit — Tr = 1 it re- 
produces the expression obtained by Blanter and Sukho- 
rukov by means of the minimal-correlation principle and 
the exact quantum- mechanical results. 

IV. CASCADE CORRECTIONS 

A straightforward extension of the minimal correlation 
approach to higher cumulants has led _to a discrepancy 
with the quantum mechanical results. tJ The reason is 
that the cavity is not just a reservoir with a noncqulib- 
rium distribution of electrons. As suggested by Eqs. (0), 
their distribution function /(e) also exhibits fluctuations. 
As the cumulants of the currents II and Ir are function- 
al of the distribution function in the cavity, its fluctua- 
tion Sf changes them too. Since the characteristic time 
scale for 5f is of the order of the dwell time of an electron 
in the cavity, these changes are slow on the scale of the 
correlation time of extraneous currents, and therefore the 
cumulants of these currents adiabatically follow 5f. This 
results in additional correlations, which may be termed 
"cascade" because lower-order correlators of extraneous 
currents contribute to higher-order cumulants of mea- 
surable quantities. One can write for the low- frequency 
transforms of the corresponding quantities 



d£ ^iw 5f{£) 



(ii) 



where S(...)/5f denotes a functional derivative of the 
corresponding quantity with respect to f(e). For exam- 
ple, the third cumulant of the current may be written as 
the sum of the minimal-correlation value 



((Il))r 



g 3 rM))-gI((i r )) 



(G l + GrY 



(12) 
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FIG. 2. The second cumulant and the two contributions to 
the third cumulant of current. The external ends correspond 
to current fluctuations at different moments of time and the 
dashed lines, to fluctuations of the distribution function. The 
full circle and triangle correspond to cumulants of extraneous 
currents and the empty triangle to the functional derivative 
of the second cumulant. 



(13) 



and the cascade correction 



The factor 3 is due to the fact that this equation in gen- 
eral includes three different currents and allows three in- 
equivalent permutations of them. 

The cascade corrections are conveniently presented in 
a diagrammatic formtj (see Figs. ^ and|3|). The rules 
for constructing these diagrams strongly differ from the 
ones known for Green's functions in quantum mechanics. 
The diagrams do not present an expansion in any small 
parameter and their number is strictly limited for a cu- 
mulant of a given order. All diagrams present graphs, 
whose outer vertices correspond to different instances of 
current and whose inner vertices correspond either to cu- 
mulants of extraneous currents or their functional deriva- 
tives. The number of arrows outgoing from an inner ver- 
tex corresponds to the order of the cumulant and the 
number of incoming arrows corresponds to the order of a 
functional derivative. Since the nth cumulant presents a 
polynomial of the distribution function of degree n, the 
number of incoming arrows at any inner vertex cannot 
exceed the number of outgoing arrows. Apparently, the 
difference between the total order of cumulants involved 
and the total number of functional differentiations should 
be equal to the order of the cumulant being calculated. 
As there should be no back-action of higher cumulants 
on lower cumulants, all diagrams are singly connected. 
Therefore any diagram for the nth cumulant of the cur- 
rent may be obtained from a diagram of order m < n by 
combining it with a diagram of order n — m + 1, i.e. by 
inserting one of its outer vertices into one of the inner 
vertices of the latter. Hence the most convenient way to 
draw diagrams for a cumulant of a given order is to start 
with diagrams of lower order and to consider all their 
inequivalent combinations that give diagrams of the de- 
sired order. The analytical expressions corresponding to 
each diagram contain numerical prefactors equal to the 
numbers of inequivalent permutations of the outer ver- 
tices. 

Unlike the case of a diffusive conductor, the third and 
fourth cumulants include now all possible diagrams and 
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not only those that are constructed of second-order cu- 
mulants (Fig. |J diagram a). The third cumulant is 

presented by diagrams b and c in Fig. 0. Diagram b 

presents the minimal-correlation value (|l 3 ) and diagram 
c presents the only possible cascade correction (13) ob- 
tained by combining two second cumulants. 

We are now in position to evaluate the diagrams. The 
functional derivative is easily obtained by differentiating 
Eq. (H) and substituting it into (O), which gives 



s((ij)) 



GlGr 



1 {G L [l-2f R + 2T R (f R -f)] 



(G L + Gj 
G R [l-2f L +2T L (f L -f)]}. 



(14) 



To calculate the fluctuation 5f, one has to write down 
equations (Q) and the current-conservation law in the 
energy-resolved form. This immediately gives 



Sf(e) 



Gl + Gr 



(15) 



where (Il,r)s ar e energy-resolved extraneous currents. 
Since fluctuations at different energies are completely in- 
dependent, one easily obtains that 



(Sf(s)5I L ) 



{Gl + G R ) 2 



G L ((I 2 )) £ -G R ((I 2 )) £ 



(16) 



Hence the total third cumulant, which is the sum of (|T^ 
and (113), is of the form 



«/£)> 



e 2 I 



(G L + G R f 



{(Gl+G r )[(G l + G r ) 2 



x(G 3 l + G R ) - 3(G L + G R ){T L G R + T R G 4 L ) 

+2T 2 L G R + 2T 2 R G L ] - 3G L G R [G 2 L (1- T R ) - G 2 R (1 - T L )] 
x[G 2 L (l-2T R )-G 2 R (l-2T L )}}. (17) 

In the case of perfectly transparent leads = T R ; = 1 
the cascade correction to the third cumulant is zero and 
the minimal-correlation result 

//t3\\_ _ 2j GlG R {Gl ~ G R ) 2 

{G L + G R Y 

is reproduced. This is why the discrepancy between the 
minimal-correlation and quantum-mechanical results was 
noted by Blanter and co-workers only for the fourth cu- 
mulant. 

The fourth cumulant is presented by a sum of six dia- 
grams shown in Fig. |3|. Diagram a presents the minimal- 
correlation value 
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t 
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FIG. 3. The contributions to the fourth cumulant of the 
current. Dashed lines correspond to fluctuations of the distri- 
bution function. Full circles, triangles and squares correspond 
to the second, third, and fourth cumulants of extraneous cur- 
rents. The empty triangles and squares present their func- 
tional derivatives. 



((Il))rn = 



G%((lt)) + Gi((I R )) 



(G L 



g r Y 



(18) 



The rest of diagrams are obtained by combining the di- 
agrams for the third and second cumulants. Diagrams 
b and c are obtained by inserting the second cumulant 
into the diagram c for the third cumulant, and diagram 
d is obtained by plugging diagram c for the third cumu- 
lant into the second cumulant. Diagrams e and / are 
obtained by inserting diagram b into diagram a and di- 
agram a into diagram b. The corresponding analytical 
expressions contain numerical prefactors 1, 6, 12, 3, 6, 
and 4 that present the numbers of inequivalent permuta- 
tions of four currents entering into the cumulant. 
The first correction is given by an expression 



Ai«l!» =6 / dsi / de 2 



5 2 ({ID) 
Sf(£i)Sf(e 2 ) 

x {5f{e 1 )5I L ){5f{e 2 )5I L ). 



(19) 



The second functional derivative 



Sf( £l )Sf(e 2 ) 



-25{ £l 



GlGr(GrTl + GlTr) 



[G l + G R ) 2 



(20) 



is obtained by differentiating (^) twice with respect to 
/(e), and the two correlators in ( [l9| ) are given by (|l6|). 

The second cascade correction is given by an expres- 
sion 



de- 



SjSfje^Sh 
Sf(e 2 ) 



4 



x (5f(e 2 )SI L ), (21) 
where the first functional derivative is given by ([OJ) , 



S(Sf{ ei )6I L ) 



26(e 1 -e 2 )- 



Gj,G 



L^R 



Sf(e 2 ) yi ZJ {G L + G R f 

+ (i-r L )f L -(i-r R )f R ], 



[(r £ - r R )f 



(22) 



and the last correlator is given by (|16|). 

The third contribution is given by an expression 



A 3 (U1)} =3 de x de 



OJ{£l) 0f{£ 2 ) 

(23) 



where the functional derivatives are given by ( |l4| ) and 
the second cumulant of the distribution function 

{ 5f{e 1 )8f{e 2 )) = e 2 5(e 1 -e 2 ) { -^%^^ (24) 

(G L + G R ) Z 

is obtained by multiplying equations ( |l5| ) with s = E\ 
and e — e 2 . 

The fourth cascade correction involves third-order cu- 
mulants of extraneous currents and is given by an expres- 
sion 

A 4 «T|» = 6 / de S 4ff^(Sf(e)6ll) m , (25) 



where 



(5f(e)Sll) m = ^ e 9MlMl±pmk (26 ) 

is obtained by multiplying one equation ( |l5| ) and two 
equations (0) and averaging them with the correlators 



The fifth correction is given by 

Wl))m 

Sf(e) 



A 5 ((li))=4j de- 



(5f(e)8I L ). (27) 



The functional derivative in the integrand 
5((I a L )) m G L G 



L Lll m _ ^L^R 

TfW~ e {G L + G R f 



{g|[1- 61^/(1-/) 



-er L (i-r L )(/-/ L ) 2 ] 
- G\ [i - 6r fl /(i - /) - 6r fl (i - r R )(f - f L ) 2 ] } 

(28) 

is calculated similarly to (|lj). 

The total fourth cumulant of current is given by the 
sum of its minimal-correlation value (nSyi and the cascade 
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FIG. 4. The circuit theory representation of the junction, 
with matrix Greens functions g~L,gR and g, matrix currents 
Ir and II and the counting field \ shown. 

corrections given by @, @, ©, @, and @. The 
full resulting expression of rather complicated form is 
given in the Appendix, and here we give only its limiting 
values 



r<2 r<2 

U L^R 



(G L + G R ) 
— 8GlG r + Cxjj] 
in the high-transparency limit Tl = F R — 1 and 



7 \pL - &GrG l + V2G L G R 

(29) 



^-" V iC&GW [Gl -" GlGR 
+ 31G|G^ - 40G 3 L G 3 R + 31G L G R - 8G L G R + G R ] 

(30) 

in the low-transparency limit where — > and T R — > 0. 

V. CIRCUIT-THEORY RESULTS 

We now show that the same results for the cumulants 
can be obtained within an ensemble-averaged quantum 
mechanical Green's function approach. We use the cir- 
cuit theory of full counting. statistics, recently developed 
by Nazarov and Bagrets,E3 which allows us to calculate 
cumulant by cumulant in a systematic way. For short- 
ness of the presentation, we only calculate the first three 
cumulants. 

The circuit theory representation of the junction is 
shown in Fig |J. It consists of three "nodes" , the two 
reservoirs and the dot, connected via two "resistances", 
the left and right point contact. Each node is assigned 
a 2 x 2 matrix Green's function, i.e. 9L,gR and g. 
The matrices are subjected to a normalization condition 
9l = 9% = S 2 = 1- The Green's functions of the two 
reservoirs are known, 



9L 



9l,r 



e lxf * /2 g L e 



1 - 
-2(1 



Vl,r(£) 



9r — 9r, 

-2/i,Jt(e) 
h,R(e)) 2/£, fl (e)-l 



(31) 



where f z is the Pauli matrix, /^^(e) are the Fermi dis- 
tribution function of the reservoirs and x is t ne counting 
field (due to current conservation, it is only necessary to 
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count the electrons in one reservoir). The matrix g(x) is 
determined from a matrix current conservation equation 



0, I 



L,R 



GL,R[9L,R,g] 

A + T L , R [{g L ,g}-2] 



(32) 



where [..,..] is the commutator and {■■,■■} the anti- 
commutator. Knowing g(x), the full counting statistics 
of charge transfer can be found. However, in the system 
under study, it is not possible to find an explicit expres- 
sion for g(x) in the general case (arbitrary N^,N R and 
Fl, Tfi)j and the full counting statistics has to be studied 
by numerical means. Here we are only interested in the 
first three cumulants, which can be found analytically 
by an expansion of the Green's functions in the counting 
field x- 

The first three cumulants are given by (evaluated at 
the left contact) 



dE tr [f z I L ] 



{{ID) = ~i 



{{iD) = 



dE tr 



dE tr 



dlL 
dx . 
d 2 I L 
dx 2 



(33) 



To calculate these cumulants, we thus need to expand 
the Green's functions, and correspondingly the matrix 
currents, to second order in the counting field x, i-e. 



six) 



X_ 
2 



j(2) S (n) 



d^l 

dx n 



(34) 



and similarly for the other quantities. For simplicity, 
we consider the case with zero temperature and the left 
reservoir held at a finite voltage eV. In this case, only 
the energies < e < eV need to be considered, where 
= 1 and /r(e) = 0, and we drop the energy nota- 
tion below. 

For the first cumulant, we have the matrix currents to 
zeroth order in the counting field, 



f(0) _ 

X L,R ~ 



4 + r L , fl [{5i° ) ii ,5 (0) }-2]' 



(35) 



where from Eq. (|3l|), one has g^ RL = g R ,L- From the 



matrix current equation in Eq. 
we then obtain 



1-2/ -2/ 
-2(1-/) 2/-1 



r(0) 



G L 



If =0, 



G L + G 



R 



(36) 



where / is the distribution function in the dot, as in Eq. 
(|l|). Knowing g(°) we find from Eq. fl35| ) and then 
the current from Eq. ( |3^ ) 

For the second cumulant, we need to expand the matrix 
currents to first order in x- Noting that the expressions 



in the matrix denominator appearing in the expansion, 
is 4 + T L<R [{gW R , g(°)} - 2] = 4, we obtain 

I^ = ^([^\9^] + [^\-9^}) 

+ §({^ 1) ^ (0) } + {^ ) ,5 (1) }), 
.£ } = %$^ W ] + ^{$\s< l) }- (37) 
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In addition to the matrix current equation 1^ + I R ^ = 0, 
we get an extra condition for g^ from the normalization 
condition g(x) 2 = lj namely 

fix) = 1 + X{.9 (0) ,5 (1) } + O(x 2 )-1 

^{.9 (0) ,.9 (1) } = (38) 

Staring from the ansatz (equivalent to the parametriza- 
tion in Ref. O) 



h 



(i) 



n 

h {1) -h 

"21 "ll 



12 
(1) 



(39) 



Eq. (H) gives h£> = - 2/)// - - /)//. In- 



serting g^ into the matrix current equation 1^ + I R 
gives h$ = + 4/ 2 and then 



h (1) 

"12 



AG, 



[G L 
GlG 2 , 



GrY 

-Gld-T, 



[Gl + G 2 L G R (1 + T R ) 



(40) 



From h\J we thus obtain all components of g^ . Insert- 
ing this into the matrix currents in Eq. ( p7| ) we get the 
second cumulant from Eq. (|33j). It coincides exactly with 
Eq. ©. 

The calculation of the third cumulant proceeds along 
the same lines. One first expands the matrix currents 
to second order in x (not presented due to the lengthy 
expressions). The requirement that the 0(% 2 ) term in 
Eq. (|f) should be zero gives {s (2) ,s (0) } + 2(#W) 2 = 0. 
Using the ansatz 



9 {2) = 



h {2) h {2) 

"ll "12 
"21 "ll 



(41) 



one gets h$ 



W\2 



(9 ( a>) 



/f + h<$(l 



h\ 2 ; (1 — /)//. The matrix current equation I L 



2f)/f 
F(2 



= 



(2) 

then gives h\-{ 



(2) 

and h 12 (not written out), which fully 
determines g^ . Inserting this into the expression for 
the matrix currents we find the third cumulant from Eq. 
d33| ) . It coincides exactly with Eq. ( |l7| ) . 

We point out that it is in principle possible to ob- 
tain analytical expressions for all higher cumulants in the 



G 



same way, although the procedure is rather cumbersome 
already for the third cumulant. 

The third and fourth cumulants of current miw be 
also obtained by means of random-matrix theoryEj us- 
ing the diagrammatic technique proposed by Brouwer 
and Beenakker J9 Substituting the resulting transmission 
probabilities for the whole system and using Eq. (||), 
one obtains expressions that coincide with Eq. ( |l7|) and 
the expression for the fourth cumulant given in the Ap- 
pendix. 



-(92r L - 60r| - 30 - 168r L r fl + 120r 2 L T R + 55T R )G R G 3 L 
+ (45r L + 30r| - 10 - 60T 2 L )G 7 R G 2 L 
+ (-36r| - 30r L + 5 + 60r|)G|Gi 

+(-i + r L )(er| - er L + i)G 9 R ] Ug l + g r ) 9 . 



VI. CONCLUSION 



In summary, we have shown that the semiclassical cas- 
cade approach gives the same results for the third and 
fourth cumulants of current in a chaotic cavity with im- 
perfect leads as the circuit theory. This leads us to the 
conclusion that this approach may be applied to a wide 
class of systems that may include both semiclassical and 
quantum-mechanical elements. The advantage of the cas- 
cade approach is its physical transparency and a relative 
simplicity. For example, if the system consists of a num- 
ber of cavities connected by contacts whose cumulants 
of current are known, this approach allows one to easily 
construct the cumulants of the current for the whole sys- 
tem. In principle, it also allows an inclusion of inelastic 
scattering processes and a calculation of cross-correlated 
cumulants of current in multiterminal systems. Therefore 
it presents a reasonable alternative to the full counting 
statistics based on the circuit theory. 
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APPENDIX 



In the case of arbitrary transmissions of the contacts 
the fourth cumulant is given by the expression 



((It)) = -e 3 / (T R ~ \)(%Y\ - 6r fl + \)G\ 



+(eor^ - 30r fl - Z%Y\ + 5)G R G 8 L 
+ (S0T% - 10 - 60rfj + A5T R )G 2 R G 7 L 

-(i2or L r^-6or^-30+92r J? +55r L -i68r L r fl )G 3 ^G ! | 

h(72r L r fl + 4 + 72T 2 R - 6T L - 96T L T 2 R - 51T R )G%G 5 L 
K4 - 6T R - 5ir L + 72r| + 72T L T R - 96T 2 L T R )G R G 4 L 
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